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ARTICLE INFO ABSTRACT 


Keywords: Porcine enteric alphacoronavirus (PEAV) is a newly identified swine enteropathogenic coronavirus that causes 
Porcine enteric alphacoronavirus watery diarrhea in newborn piglets. In this study, an original, highly virulent PEAV strain GDS04 was serially 
Attenuation passaged in Vero cells. The virus titers and sizes of syncytia increased gradually with the cell passages. Newborn 
Pathogenicity 


piglets were orally inoculated with PEAV P15, P67 and P100. Compared with P15 and P67, P100 resulted in only 
mild clinical signs and intestinal lesions in piglets. The virus shedding in feces and viral antigens in intestinal 
tract were markedly reduced in P100-inoculated piglets. Importantly, all P100-inoculated newborn piglets 
survived, indicating that P100 was an attenuated variant. Sequence analysis revealed that the virulent strain 
GDS04 had four, one, six and eleven amino acid differences in membrane, nucleocapsid, spike and ORFlab 
proteins, respectively, from P100. Furthermore, more differences in the predicted three-dimensional structure of 
S protein between GDS04 and P100 were observed, indicating that these differences might be associated with the 
pathogenicity of PEAV. Collectively, our research successfully prepared a PEAV attenuated variant which might 


Newborn piglets 
Genomic analysis 


serve as a live attenuated vaccine candidate against PEAV infection. 


1. Introduction 


Porcine enteric alphacoronavirus (PEAV) was first detected by our 
team via genomic analysis of samples collected from a diarrhea-out- 
break swine herds routinely vaccinated with porcine epidemic diarrhea 
virus (PEDV) vaccine in a farm in Guangdong, China in February 2017 
(Gong et al., 2017). This novel swine enteric coronavirus emerged in 
China at least since August 2016 and it has been widely detected in 
Guangdong, China, with a prevalence of 43.53% by a retrospective 
detection study (Zhou et al., 2019). Recently a PEAV-like strain CN/ 
FJWT/2018 was recently discovered in Fujian, China (Li et al., 2018), 
indicating a prevailing trend in pig farms. 

Since the first report, the complete genome of the PEAV strain 
GDS04 was sequenced (Gong et al., 2017). PEAV is an enveloped, 
single-stranded, positive-sense RNA virus with a genome of appro- 
priately 27 kb in length (Gong et al., 2017). The genome organization of 
PEAV is similar to that of bat-like HKU2 strains of coronavirus, with an 
order of: 5’ untranslated region (UTR), open reading frame 1a/1b 
(ORF1a/1b), spike (S), nonstructural protein 3 (NS3), envelope (E), 


membrane (M), nucleocapsid (N), nonstructural protein 7a (NS7a), and 
3’ UTR (Lau et al., 2007). The S protein of coronaviruses (CoVs), the 
pivotal surface glycoprotein, is involved in virus attachment and entry, 
and induction of neutralizing antibodies in vivo (Cruz et al., 2008; Woo 
et al., 2010). GDS04 strain of PEAV has the smallest S protein among all 
CoVs (Gong et al., 2017). 

The clinical symptoms caused by PEAV in newborn piglets are si- 
milar to that by other porcine enteric pathogens such as PEDV and 
transmissible gastroenteritis virus (TGEV), which include vomiting, 
diarrhea, dehydration, and mortality rate as high as 90% (Gong et al., 
2017; Zhou et al., 2018b). Since PEAV was detected (Gong et al., 2017), 
another two swine enteric HKU2-related CoVs named SADS-CoV and 
SeACoV were identified in the same region, causing same diarrheal 
disease as PEAV strain GDS04 by experimental infection (Pan et al., 
2017; Xu et al., 2019a; Zhou et al., 2018b). 

Currently, live-attenuated vaccine has been widely used to prevent 
and control porcine enteric pathogens, such as PEDV, with good clinical 
effect (Song et al., 2007). Live-attenuated vaccine is prepared from less 
virulent strain. The virulence of pathogenic microorganisms was 
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reduced by various methods (Blanco-Lobo et al., 2019; Jie et al., 2018; 
O’Donnell et al., 2015), such as genetic manipulation (O’Donnell et al., 
2015), and cell passage is the most commonly used (Jie et al., 2018). 
Usually, the virulent parent strain was continuously inoculated cell 
lines more than 100 generations in vitro, and in this process, some 
mutations will happen in the viral genome, which might lead to the 
strain's virulence reduce. These cell-passaged strains are then tested in 
animals to determine whether they have less virulence and higher im- 
munogenicity before being used to prepare live-attenuated vaccines. Of 
note, several groups have used cell passage to prepare porcine CoVs 
attenuated strains, such as PEDV (Lin et al., 2017) and TGEV (Motovski 
et al., 1985), which laid a good foundation for preparation of live-at- 
tenuated vaccines. 

PEAV is an important enteropathogen in pigs, but currently no ef- 
fective treatments or vaccines against PEAV infection are available. To 
develop a live-attenuated vaccine for PEAV, we generated candidates 
via serial cell passage of the parental PEAV GDS04 strain and evaluated 
the pathogenicity of PEAV P15, P67 and P100 in 5-day-old newborn 
piglets. Further, we identified genetic changes related to attenuation by 
performing comparative, complete genomic sequence analysis of viru- 
lent and attenuated PEAV variants. 


2. Materials and methods 
2.1. Serial virus passaging in Vero cells 


Vero cells were obtained from ATCC (ATCC number: CCL-81) 
(USA), and cultured in Dulbecco' s modified eagle medium (DMEM) 
(Hyclone, USA) supplemented with 100 U/mL penicillin, 100 U/mL 
streptomycin, and 10% fetal bovine serum (FBS) (BOVOGEN, 
Australia). The maintenance medium for PEAV propagation was DMEM 
supplemented with 10 g/mL trypsin (Gibco, USA). The isolation and 
identification of PEAV GDS04 strain were reported previously by our 
laboratory (Xu et al., 2019a). 

Serial virus passaging and plaque-purification in Vero cells were 
performed as previously described with some modifications (Lin et al., 
2017; Xu et al., 2019a). Briefly, Vero cells were cultured in 6-well 
plates, and washed three times with sterile pH 7.4 1 x phosphate buffer 
saline (PBS) when 90% confluent. Ten uL of PEAV GDS04 strain to- 
gether with 2mL maintenance medium were added to plate and cul- 
tured continuously at 37°C in 5% CO, to observe cytopathic effect 
(CPE). Around 1 day post inoculation (d.p.i.), the plates were frozen at 
-80°C and thawed twice after CPE was evident in the inoculated cell 
monolayers. The cells and supernatants were harvested together and 
used as seed stocks for the next passage. PEAV GDS04 strain was pas- 
saged regularly to a total of 100 passages in Vero cells every 2-3 days 
and the P15, P67 and P100 were plaque-purified. Finally, T175 flasks 
were used for the P15, P67 and P100 propagation in Vero cells. 


2.2. Pigs 


Four-day-old crossbred conventional newborn piglets produced by 
Duroc x Landrace x Native pigs of Guangdong of China were procured 
from Wen' s Foodstuffs Group Co, Ltd (Guangdong, China). All pigs 
were maintained in our animal facility with a mixture of skim milk 
powder (Inner Mongolia Yi Li Industrial Group Co., Ltd, China) with 
warm water ad libitum for 1 day before the experimentation. The an- 
imal study was approved by the Institutional Animal Care and Use 
Committee of the Sun Yat-sen University (Guangdong, China) and an- 
imals were treated in accordance with the regulations and guidelines of 
this committee. 


2.3. Measurement of virus growth in Vero cells 


Kinetics of PEAV replication in Vero cells were performed as pre- 
viously described with some modifications (Lin et al., 2017). Briefly, 
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Vero cell monolayer was inoculated with P15, P67 and P100 at a 
MOI = 0.1. The virus inoculum was removed at 2h post inoculation 
(hpi), washed with pH 7.4 1 x PBS twice and replaced by the main- 
tenance medium. Culture supernatant and cell lysates were collected at 
6, 12, 18, 24 and 30 hpi, respectively. After freezing and thawing once, 
the mixtures of culture supernatants and cell lysates were subjected to 
titration for 50% tissue culture infectious dose (TCIDs9) in 96-well 
plates and calculated using Reed-Muench method (Reed and Muench, 
1938). 


2.4. Experimental infection of PEAV at selected passages in conventional 
newborn piglets 


Thirty-two 4-day-old conventional newborn piglets were randomly 
divided into four groups and 8 piglets in each group. Prior to inocula- 
tion, piglets were confirmed to be negative for the major porcine enteric 
viruses including PDCoV, PEDV, TGEV, PRoV, PEAV, by testing of 
rectal swabs using specific RT-PCR according to previously described 
method (Xu et al., 2018; Zhou et al., 2018a). In addition, serum of each 
piglet was collected and a commercial kit (ALLRIGHT BIOTECHONO- 
LOGY, China) was applied in order to detect PEDV-specific antibody 
according to the manufacturer' instruction. PEDV-specific antibodies 
were not detected in all sera samples (data not shown). After 1-day 
acclimation, piglets in group 1 were orally inoculated with 5mL of 
maintenance medium and served as uninfected controls. Piglets in 
group 2, 3 and 4 were orally challenged with 5mL of maintenance 
medium containing a total of 5 x 10° TCIDso of P15, P67 and P100, 
respectively. To reduce the risk of cross contamination among piglets 
and piglet death caused by non-experimental factors, the entire animal 
study period was ended in nine days after challenge. 

All piglets were observed daily for clinical signs of diarrhea and 
lethargy. Diarrhea severity was scored with the following criteria (Chen 
etal., 2015): 0 = normal, 1 = soft (cowpie), 2 = liquid with some solid 
content, 3 = watery with no solid content. 

Rectal swabs were collected daily from each piglet from 1 d.p.i. to 9 
d.p.i., submerged into 1 mL sterile pH 7.4 1 x PBS immediately after 
collection and stored at -80°C for virus shedding analysis. Two piglets 
from each group were necropsied at 7 d.p.i.. At necropsy, the fresh 
jejunums were collected and then formalin-fixed. The formalin-fixed 
samples were used for histopathology and immunohischemistry ana- 
lysis. In addition, the mortality of newborn piglets in each group was 
recorded daily. 


2.5. RNA isolation and real-time RT-PCR analysis 


Viral RNA shedding in piglet feces after infection with PEAV was 
detected as previously described with some modifications (Xu et al., 
2019a). Briefly, total RNA was prepared from the supernatant of rectal 
swab from each piglet using a RNeasy kit (Magen, China) according to 
the manufacturer’ instruction, and was treated with DNase I. Two ug of 
total RNA was used for cDNA synthesis by reverse transcription using 
RT-PCR kit (TaKaRa, Dalian). The specific primers for the n gene of 
PEAV (sense: 5’-CTGACTGTTGTTGAGGTTAC-3’; antisense: 5’-TCTGCC 
AAAGCTTGTTTAAC-3’), and probe (5’-FAM-TCACAGTCTCGTTCTCGC 
AATCA-TARMA-3’) were designed with reference to the previous pub- 
lication (Zhou et al., 2018a) and synthesized by Invitrogen Company 
(Shanghai, China). The real-time PCR assay was carried out with an 
Applied Biosystem 7500 instrument (Life Technologies, USA). The PCR 
was performed in a 20-uL volume containing 1 uL of cDNA, 10uL of 
Thunderbird Probe qPCR Mix, 0.04uL 50 x Rox reference dye 
(TOYOBO, Shanghai), 0.2 umol/L of probe, and a 0.3 mol/L of each 
gene-specific primer. The thermal cycling parameters were as follows: 
95°C for 30s; 45 cycles of 95°C for 5s, 62°C for 30s. And the n gene 
was amplified to generate the standard curve from PEAV GDS04 strain 
using the specific primers (sense: 5’— CCGCTCGAGATGGCAACTGTTA 
ATTGG-3’; antisense: 5’-CGCGGATCCCGATTAATAATCTCATCCAC-3’) 
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Fig. 1. The flowchart of the isolation and passaging of PEAV GDS04 strain in Vero cells. 


that were designed with reference to the published sequence (GenBank, 
Accession no: MF167434.1), and the PCR products were cloned into the 
pEGFP-N1 (Clontech, USA). The known plasmid concentration was 10- 
fold serially diluted for generating a standard curve in each plate. The 
quantity of PEAV viral RNA in tested samples was calculated based on 
the cycle threshold (Ct) values for the standard curve. 


2.6. Histological and immunohistochemical staining 


Histological and Immunohistochemical staining were performed as 
previously described with some modifications (Xu et al., 2019a). 
Briefly, tissue samples of jejunum of the piglets from the challenged and 
control groups were separated, routinely fixed in 10% formalin for 36 h 
at room temperature, dehydrated in graded ethanol, embedded in 
paraffin, cut in 5-um sectioned, and mounted onto glass slides. After the 
sections were deparaffinized, rehydrated, and stained with hematoxylin 
and eosin (H&E), the slides were examined and analyzed with con- 
ventional light microscopy. Five-um sections of formalin-fixed paraffin- 
embedded tissues were placed onto positively charged glass slides and 
the slides were air dried for 120 min at 60°C. The tissue sections were 
deparaffinized, rinsed and incubated with target retrieval solution 
(Servicebio, China). After being blocked with 1% BSA (Solarbio, China), 
the sections were incubated with anti-PEAV M monoclonal antibody 
(Wen's Foodstuffs Group Co., Ltd, China) as the primary antibody and 
the dilution of antibody was 1:200 for 12h at 4°C. They were then 
incubated with peroxidase-labeled goat anti-mouse IgG secondary an- 
tibody (Dako, Denmark) and the dilution of antibody was 1:200 for 
50 min at room temperature, and the samples were finally visualized 
with a 3, 3’-diaminobenzidine (DAB) chromogen kit (Dako, Denmark). 
Hematoxylin was used for counterstaining. Tissues of piglets from un- 
infected control groups were used as negative samples. The im- 
munohistochemistry slides were evaluated according to the evaluation 
system of histology and immunohistochemistry by Jung et al. (2014). 


2.7. Genomic sequencing and sequence analysis 


The complete genomes of P15, P67 and P100 were sequenced by 
IlluminaHiSeq as previously described with some modifications 
(Malboeuf et al., 2013). Briefly, total RNA was prepared from the PEAV- 
infected Vero cells using a RNeasy kit (Magen, China) according to the 
manufacturer’ instruction, and was purified as described above. cDNA 
synthesis was performed by reverse transcription using RT-PCR kit 
(TaKaRa, Dalian). The genomic cDNA was disrupted and converted to a 
sticky end by adding the A-base at 3’ end of the cDNA. The DNA con- 
taining Index sequence was added at either side of sticky end by 


complementary bases. The target fragment within the scope of a certain 
length was collected before magnetic beads selection. The sequencing 
library was built, tested using PCR amplification Index sequence, and 
combined to the chip by bridge PCR. At last, the complete genome was 
sequenced by Ilumina HiSeq. A genome homology analysis and phy- 
logenetic trees were constructed by using the maximum likelihood 
method with MEGA 5 software (http://www.megasoftware.net/) based 
on the whole-genome nucleotide sequences from PEAV GDS04, P15, 
P67 and P100 together with other CoVs, like PEDV, PDCoV and TGEV. 
Sequence alignments of N, M, S and ORFlab amino acids (aa) of GDS04, 
P15, P67 and P100 were performed using the DNAMAN. 


2.8. Homology modeling of the S proteins of P15, P67 and P100 


The S protein of CoV, the pivotal surface glycoprotein, plays an 
important role in virus attachment and entry, and induction of neu- 
tralizing antibodies in vivo (Cruz et al., 2008; Li et al., 2017, 2016). To 
further determine the differences of Vero-passaged PEAV strains, the 
three-dimensional structures of S proteins of GDS04, P15, P67 and P100 
were predicted in SWISS-MODEL using amino acid homology modeling, 
based on the existing three-dimensional structure of S glycoprotein 
trimer of Human coronavirus NL63 (PDB code 5SZS) (Walls et al., 
2016) in Protein Data Bank. The PyMOL Molecular Graphics System 
(DeLano Scientific; http://www.pymol.org) was used for figures pre- 
paration. 


2.9. Statistical analysis 


Statistical comparisons were performed using GraphPad Prism 
software. The significance of the differences between the growth rate of 
P100 and P15 or P67 in the TCIDsy was determined by the ANOVA and 
Mann-Whitney accordingly. 


3. Results 
3.1. PEAV P100 had highest infectious titers in Vero cells 


To generate an attenuated PEAV vaccine candidate, PEAV GDS04 
was passaged regularly to a total of 100 passages in Vero cells and the 
selected passages including P15, P67 and P100 were characterized by 
sequencing and analyzing the complete genome (Fig. 1). To determine 
the growth rate of P15, P67 and P100 in vitro, we inoculated P15, P67 
and P100 at MOI = 0.1 in Vero cells and detected the live virus count at 
corresponding time point with a TCIDso assay. Interestingly, we ob- 
served that the sizes of PEAV-induced syncytia increased in P67 and 
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Fig. 2. The growth kinetics of PEAV P15, P67 
and P100 in Vero cells. 

(A) Mock-inoculated Vero cell culture re- 
presenting normal cells. (B-D) Vero cells were 
inoculated with PEAV P15, P67 and P100 at 
MOI = 0.1. PEAV-inoculated Vero cells at 
24 hpi. showing syncytium and cells detach- 
ment and was indicated by arrows. (E) Cell 
lysates were sampled at designated time points 
and titrated using TCIDso infectivity assay. 
Results are representative of three independent 
experiments. Data are represented as 
mean + SD, n = 3. *stands for p < 0.05. 
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P100 (Fig. 2A). The diameter of P100-syncytia was larger than that of 
P15 and P67 (Fig. 2A). In addition, the muti-step growth kinetics of 
PEAV in Vero cells showed similar growth curve trends of P15, P67 and 
P100, but the titers of P100 were significantly higher than those of P15 
and P67 at 24h.p.i. (P < 0.05) (Fig. 2B). 


3.2. PEAV P100 demonstrated reduced clinical signs and fecal virus in 
newborn piglets 


To evaluate the pathogenicity of Vero-passaged PEAV strains in 
newborn piglets, we experimentally infected 5-day-old conventional 
newborn piglets with P15, P67 and P100 at a dose of 5 x 10° TCIDso/ 
head via oral feeding. As shown in Fig. 3A, the newborn piglets in- 
oculated with P15, P67 or P100 all showed diarrhea from 1 d.p.i. to 9 
d.p.i., as compared with the control. But compared to P15 and P67, the 
degree of diarrhea of the newborn piglets inoculated with P100 was the 
lowest, indicating that viral pathogenicity decreased gradually in 5- 
day-old newborn piglets with cell passage. In addition, we explored the 
fecal viral shedding in PEAV-inoculated piglets from 1 d.p.i. to 9 d.p.i. 
by qRT-PCR. Consistent with the clinical signs, the amount of virus RNA 
was lower in P100-inoculated newborn piglets, and no PEAV RNA was 
detected in the negative control piglets during the study (Fig. 3B). 
Taken together, these results demonstrated that cell passage gradually 
decreased viral pathogenicity of GDS04 to 5-day-old newborn piglets. 


3.3. PEAV P100 was attenuated in newborn piglets 


Since PEAV at selected passages caused varying degrees of diarrhea 
in newborn piglets, we also recorded the mortality of newborn piglets 
in each group. As shown in Fig. 4, in P15-inoculated group, except two 
piglets that were necropsied at 7 d.p.i., 2, 1, 1, and 1 death (s) were 
recorded at 3, 4, 5, and 8 d.p.i., respectively. Similarly, except two 
necropsied piglets, 4 piglets in P67-inoculated group died from 5 to 9 
d.p.i. No piglets died in P100-inoculated group and the control group 
except two piglets that were necropsied at 7 d.p.i.. Taken together, 
these results suggest that PEAV P100 strain is low pathogenic to the 
newborn piglets. 


3.4. PEAV P100 caused reduced lesions and had limited antigen distribution 
in intestinal tract 


To determine the gross pathological and histological changes in 
piglets infected with PEAV GDS04 at selected passages, two piglets from 
each group were necropsied at 7 d.p.i. Gross findings were similar in 
piglets orally inoculated with P15 and P67. The whole intestinal tract, 
where yellow watery contents accumulated, were transparent, thin- 
walled, and gas-distended (Fig. 5B&C). No lesions were observed in the 
whole intestinal tract of the P100-inoculated piglets and the negative 
control piglets (Fig. 5A&D). Microscopic lesions were also analysed. As 
shown in Fig. 5F&G, abruption of intestinal villus was observed in P15 
or P67-inoculated piglets, whereas the intestines in P100-inoculated 
piglets and the negative control were normal (Fig. 5E&H). Consistent 
with the histopathological results, PEAV antigen was detected in the 
cytoplasm of the villous enterocytes of the P15 or P67-challenged 
piglets by immunohistochemical analysis, but no antigen observed in 
P100-inoculated and mock piglets (Fig. 5I-L). Taken together, these 
results indicate that PEAV P100 is an attenuated variant. 


3.5. Phylogenetic analysis of complete-genome and comparison of the aa 
changes among M, N, S and ORF1ab proteins 


To determine the cause of reduced pathogenicity of Vero-passaged 
PEAV strains, the complete genomes of P15, P67 and P100 were se- 
quenced and analyzed. Compared to the complete genome of PEAV 
GDS04 (MF167434.1), P15 possesses 3 point mutations, a 3-nt deletion 
and a 3-nt insertion; P67 possesses 18 point mutations, a 3-nt deletion 
and a 3-nt insertion; P100 possesses 18 point mutations, a 3-nt deletion, 
a 6-nt deletion and a 3-nt insertion, respectively (data not shown), 
presenting a clue of different pathogenicity of these three strains. In 
addition, the parental strain GDS04 shares 99.97%, 99.91%, 99.89% 
nucleotide identity with P15, P67 and P100, respectively (data not 
shown), indicating that differences in viral genome increased gradually 
with cell passages. In spite of some differences in viral genome of these 
strains, phylogenetic analysis showed that PEAV GDS04, P15, P67 and 
P100 were clustered into a large clade (Fig. 6A). We further analyzed 
the aa of proteins which might be associated with the pathogenicity of 
CoVs, such as M (Zhang et al., 2016), N (Hu et al., 2017) and S proteins 
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Fig. 3. Clinical signs and fecal viral shedding in newborn piglets inoculated with PEAV P15, P67 and P100 via oral feeding. 
(A) Average diarrhea scores after P15, P67 and P100 infection. (B) Ct values and copies of PEAV RNA of feces from piglets after P15, P67 and P100 or mock 


inoculation. 


(Hou et al., 2019b). As shown in Fig. 6B, compared to the M protein of 
PEAV GDS04, P15 possesses 2 aa mutations, P67 possesses 3 aa muta- 
tions, and P100 possesses 4 aa mutations. Compared to the N protein of 
PEAV GDS04, interestingly, only P100 possesses 1 aa mutation, while 
no mutations were observed in P15 and P67. Compared to the S protein 
of PEAV GDS04, P15 possesses 1 aa mutation, P67 possesses 3 aa mu- 
tations, while P100 possesses 3 aa mutations and 3 aa deletions. In 
addition, compared to the ORFlab protein of PEAV GDS04, P15 pos- 
sesses 2 aa mutations, 2 aa insertions, P67 possesses 10 aa mutations, 2 
aa insertions, P100 possesses 8 aa mutations, 1 aa deletion and 2 aa 
insertions. The detailed information of aa changes might help to explain 
the low pathogenicity of P100 to the newborn piglets. 


3.6. PEAV P100 has significant differences in the S protein structure 


Considering the obvious aa mutations and deletions in S proteins of 
Vero-passaged PEAV strains, we attempted to predict and analyze the 
three-dimensional structures of S proteins of these PEAV strains. After 
analysing the amino acid sequences, the S protein in GDS04 was found 
to have 27.68% identity with the S protein of Human CoV NL63 (Walls 
et al., 2016) The S protein structures of these strains were predicted by 
SWISS-MODEL, using the S protein of Human CoV NL63 as the template 
(Fig. 7). The predicted structure of GDS04 S protein forms a homo- 
logous trimer with an expect-value of 0.36, and each monomer is 
composed of 983 residues in the aa positions 31-1013 (Fig. 7A&B). In 
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Fig. 4. The survival rate of newborn piglets after infection with PEAV P15, P67 
and P100. 
The mortality of newborn piglets in each group was recorded from 1 to 9 d.p.i.. 


addition, the predicted three-dimensional structures of the S proteins of 
PEAV P15, P67 and P100 have similar overall structures when com- 
pared with GDS04 S protein (Fig. 7C-F). However the monomer struc- 
tural overlap of the four S proteins showed that there were three loops 
and two ß strands having significant structural differences in the main 
chain, especially between P100 and GDS04. 
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4. Discussion 


Since the first report of PEAV in pigs in early February of 2017 in 
Guangdong, China (Gong et al., 2017), this novel swine enteric CoV has 
been widely detected in Guangdong, China, with a prevalence of 
43.53% and an emergence in China as early as August 2016 by a ret- 
rospective detection study (Zhou et al., 2019). Although a few studies 
have demonstrated PEAV was highly pathogenic to newborn piglets 
(Pan et al., 2017; Xu et al., 2019a; Zhou et al., 2018b), there are no 
effective treatments or any vaccines against PEAV infection. In the 
present study, attenuated PEAV variants were generated by serial cell 
passaging and orally inoculated to newborn piglets to evaluate the at- 
tenuation in vivo. Further, we identified the molecular basis of at- 
tenuation by performing comparative, complete genomic sequence 
analysis of virulent and attenuated PEAV variants. 

Swine enteropathogenic CoVs affects all sectors of the pig industry 
and all cycles of production, which results in significant economic 
losses. Currently, vaccination remains the most effective measure 
against swine enteropathogenic CoVs. Of note, several groups have used 
live-attenuated vaccines to prevent and control PEDV and TGEV in pigs 
(Aynaud et al., 1991; Song et al., 2007). As a newly swine en- 
teropathogenic CoV, there are no reports about live-attenuated vaccine 
for PEAV. Pathogenic microorganisms serially passaged in cells (Jie 
et al., 2018) is one of the most common methods to prepare live-atte- 
nuated vaccine. Vero cells as a stable African green monkey kidney cell 
line (Rhim et al., 1969) is commonly used to attenuate CoVs like PEDV 
(Lin et al., 2017). In this study, PEAV GDS04 was passaged regularly to 
a total of 100 passages in Vero cells and the growth rate of P15, P67 and 
P100 was determined in vitro. The muti-step growth kinetics of P15, 
P67 and P100 in Vero cells were similar, but the titers of P100 were 
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Fig. 5. Intestinal changes in newborn piglets inoculated with PEAV P15, P67 and P100. 

(A) Macroscopic picture of intestine from a control piglet at 7 d.p.i. (B-D) Macroscopic pictures of intestine from P15, P67 and P100-challenged newborn piglets at 7 
d.p.i. Thin-walled intestinal tracts were indicated by arrows. (E) H&E-stained jejunum tissue section of a control piglet at 7 d.p.i. (F-H) H&E-stained jejunum tissue 
section of P15, P67 and P100-challenged piglet at 7 d.p.i. Blunt intestinal villus was indicated by arrows. (I) Immunohistochemically stained jejunum tissue section of 
a control piglet at 7 d.p.i. (J-L) Immunohistochemically stained jejunum tissue section of P15, P67 and P100-challenged piglet at 7 d.p.i.. 
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Fig. 6. Phylogenetic trees constructed on the basis of the complete genome sequences of PEAV P15, P67 and P100 or other CoVs and sequence alignments of PEAV at 


selected passages. 


(A) The dendrogram was constructed using the neighbour-joining method in the MEGA software package, version 5 (http://www.megasof tware.net). Bootstrap 
resampling with 1,000 replicates was performed, and bootstrap values are indicated for each node. Reference sequence obtained from GenBank is indicated by strain 
name. The scale bar represents 0.5 nucleotide substitutions per site. (B) Sequence alignments of the amino acids of M, N S and ORF 1ab proteins were performed using 


DNAMAN software. 


significantly higher than those of P15 and P67 at 24hpi, similar to 
observations in attenuated PEDV vaccine (Lin et al., 2017), indicating 
that P100 had better adaptation in cells in vitro. Usually, adaptation in 
cell lines coincides with attenuation of the viruses, such as African 
swine fever virus (ASFV) isolate Georgioa and PEDV isolate PC22A, in 
swine (Krug et al., 2015; Lin et al., 2017). The higher fitness of P100 in 
Vero cells might suggest the changed pathogenicity of PEAV. 

For a vaccine candidate, safety is the top concern. It was reported 
that swine enteric CoVs like PDCoV, PEDV and TGEV are age-depen- 
dent disease (Moon et al., 1975; Shibata et al., 2000; Xu et al., 2019b). 
Similarly, in another study, we found that PEAV was highly pathogenic 


to newborn piglets (Xu et al., 2019a), but did not cause vomiting and 
death in weaned piglets at a medium dose of 1 x 10° TCIDso within 7 
days (data not shown), for which reason we choose newborn piglets to 
conduct PEAV pathogenicity study. It was known that PEAV infection 
caused typical clinical symptoms characterized by watery diarrhea via 
oral feeding in newborn piglets (Xu et al., 2019a). In the study, we 
infected the 5-day-old newborn piglets with P15, P67 and P100 via oral 
feeding. Compared with P15 and P67, P100 caused slight diarrhea in 
newborn piglets, suggesting that the pathogenicity of PEAV changed 
through cell passage. For enteric viruses, viral shedding is often ac- 
companied by diarrhea, so was PEAV infection (Xu et al., 2019a). Thus, 
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Fig. 7. Structural analysis of the S proteins of PEAV GDS04, P15, P67 and P100. 
All the predicted structures of S protein were shown in cartoon style with different colors. Blue, yellow and purple stand for GDS04; green stands for P15; cyan stands 
for P67; orange stands for P100. (A&B) The predicted overall structure of S protein of PEAV GDS04 strain. (C-F) Three-dimensional structure of the monomer S 


protein in P15 (C), P67 (D), P100 (E) and GDS04 strain (F). (G) Structural overlap of the S proteins of PEAV at selected passages. 


we collected fecal swabs from the PEAV-challenged piglets and detected 
the fecal shedding by real-time PCR. Consistent with the clinical signs, 
the amount of virus RNA was lower in P100-inoculated newborn pig- 
lets. Interestingly and importantly, the presence of PEAV RNA in fecal 
samples of P100 inoculated piglets from 1 d.p.i to 9 d.p.i suggested the 
successful colonization and replication of P100 in intestinal tract, which 
is essential to induce mucosal immunity for an enteric pathogen vaccine 
candidate. It was reported that the lesions were only observed in in- 
testinal tract but not any other organs after PEAV infection and PEAV 
antigen was detected in the villous enterocytes of the PEAV-challenged 


piglets (Xu et al., 2019a). Compared to P15 and P67, P100 strain caused 
slight macroscopic and microscopic lesions in intestinal tract and no 
PEAV antigen was detected in the cytoplasm of the villous enterocytes 
of the P100-challenged newborn piglets by immunohistochemical 
analysis, indicating that the infectiveness of P100 was reduced. Like 
other swine enteric CoVs, such as PEDV (Shibata et al., 2000), PEAV 
infection caused fatal disease in newborn piglets (Xu et al., 2019a). 
Thus, death after infection is the key indicator to determine the at- 
tenuation of the viral strain. Compared to P15 and P67, no piglets died 
in P100-inoculated group and the control group during the study, 
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indicating that P100 is low pathogenic to the newborn piglets. Taken 
together, these results demonstrated that P100 strain is safe for pigs and 
effective in colonizing and replicating in intestine, suggesting its future 
application as a good PEAV vaccine candidate. 

Compared to P15 and P67 strains, the virulence of P100 was sig- 
nificantly reduced in 5-day-old newborn piglets. The S protein of CoVs 
is the pivotal surface glycoprotein involved in virus attachment and 
entry, virus attenuation and induction of neutralizing antibodies in vivo 
(Cruz et al., 2008; Lin et al., 2017; Woo et al., 2010). In addition, the 
virulence of a recombinant PEDV with inactivation of the endocytosis 
signal of the S protein was significantly reduced in piglets (Hou et al., 
2019a), indicating that the S protein might be associated with the 
virulence of CoVs. Compared to the S protein of parental strain GDS04, 
P100 possesses 2 aa mutations in aa positions 100 and 670, which were 
also found in P15 and P67, indicating that these 2 aa-mutations are 
insufficient to cause viral attenuation. Of note, Compared to P15 and 
P67, P100 possesses 3 aa deletions and 1 aa mutation in aa positions 
1001, 1002, 1003 and 1129, respectively, indicating that the absence 
and mutation of these amino acids might affect the viral attenuation. 
Usually, protein with complete structure has corresponding function, 
and the changes in the protein’s structure might affect its function. 
Overlap of the S protein structures of GDS04, P15, P67 and P100 
showed that the difference between P100 and GDS04 was most sig- 
nificant (Fig. 7), indicating a potential relationship between these dif- 
ferences of S protein and pathogenicity of PEAV. It was reported that 
the S protein is necessary but not sufficient for the virulence of PEDV 
(Wang et al., 2018). M protein plays an important role in the process of 
virus assembly by interacting with S and N proteins (de Haan et al., 
1999; Nguyen and Hogue, 1997). N protein can encapsulate the viral 
genome RNA to form the nucleocapsid (Spaan et al., 1983; Sturman 
et al., 1980). Inactivation of the viral 2’-O methyltransferase of nsp16 
could help to engineer a live attenuated PEDV vaccine (Hou et al., 
2019a). In this study, aa changes were also found in M, N and ORFlab 
proteins in P100 strain as compared to GDS04 (Fig. 6B). In addition, no 
additional mutation was observed by sequencing the viruses shed from 
the challenged piglets (data not shown). Therefore, these results suggest 
that the complete attenuation of a PEAV strain might be the result of a 
combination of multiple mutations along the genome, and can occur via 
multiple molecular mechanisms. Whether individual or a combination 
of genetic changes of PEAV strains alter viral infectivity, pathogenicity 
and replication efficiency need further examination by reverse genetics 
technology. Together, all these results confirmed that P100 generated 
in this study was low pathogenic to the newborn piglets and might 
serve as a good PEAV vaccine candidate. However, there are still sev- 
eral important questions needed to be addressed. For instance, what is 
the immunogenicity of P100 strain as a live-attenuated vaccine candi- 
date in pigs? What is the protective effect of P100 strain in pigs? Elu- 
cidation of these questions will help us to develop a safe and efficacious 
live-attenuated vaccine to control PEAV. 

In summary, our research successfully generated an attenuated 
PEAV variant P100. P100 has the highest titers of virus load, via serial 
passaging of the parental GDS04 strain. Remarkably, inoculation of 
newborn piglets with P100 by oral feeding caused the lowest level of 
clinical symptoms, including slight diarrhea, fecal viral shedding, and 
pathological changes with a survival rate of 100% in newborn piglets. 
The changes in genomic composition and structures found in P100 by 
genomic analysis might account for the underlying molecular me- 
chanisms of PEAV attenuation. 
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